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Abstract 

We study neutrino oscillations in vacuum and in matter using field theoretical 
methods and wave-packets. In particular, we calculate the neutrino propaga- 
tor in the presence of matter with constant density for the case of two flavors, 
obtaining the resonance formula. In the extreme relativistic limit, the result 
of the usual quantum mechanical treatment is recovered with interesting, but 
small modifications. 

I. INTRODUCTION 

Evidence for neutrinos having small, but finite, masses is accumulating fast. The main 
paradigm for this is the idea, first proposed by Pontecorvo [1], that if the neutrino fla- 
vors were a superposition of massive eigenstates, then the neutrinos could oscillate between 
these different flavors. There are many experiments looking for neutrino oscillations: solar 
neutrino measurements, atmospheric neutrino measurements, reactor experiments, and ac- 
celerator experiments. Several of them have claimed evidence for neutrino oscillations. The 
final breakthrough came in June 1998 when the Super-Kamiokande collaboration reported 
strong evidence for neutrino oscillations from the atmospheric neutrino UP-DOWN asymme- 
try. The measurements on the depletion of atmospheric muon neutrino flux fit well to a two 
flavor <-» v T oscillation model. This summer the SNO collaboration [2] showed that the 
solar neutrino deficit, pioneered by Davies [3], was due to conversion of electron neutrinos 
to mu- and tau neutrinos. Recently, also the KamLAND experiment [4] shows evidence for 
neutrino oscillations. 

In the standard quantum mechanical treatment of neutrino oscillations, the mass eigen- 
states are assumed to be relativistic and to have the same momentum, and thus, different 
energies. The familiar quantum mechanical model describing the flavor mixing process has 
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several conceptual difficulties compared to quantum field theory models, see Refs. [5-9]. For 
example, energy momentum conservation in the production and detection processes that 
forces neutrinos to be in a mass eigenstate is incompatible with neutrino oscillations. The 
neutrino oscillation probability is independent of the details concerning the production and 
detection processes only in the extremely relativistic limit. Hence, in the case that some 
of the mass eigenstates cannot be considered to be extremely relativistic, one has to use 
quantum field theory. Of course, the quantum field theory expression must reproduce the 
quantum mechanics oscillation probability in the ultra- relativistic limit. A very detailed 
review regarding different aspects and questions of neutrino oscillations in quantum field 
theory can be found in Ref. [5]. Another interesting question which is important in both a 
quantum mechanical or quantum field theory treatment is if there exists a Fock space for 
the flavor eigenstates, since there exists a Fock space for the mass eigenstates. This question 
has been discussed by Giunti et al. [10]. Fuji et al. [11] give arguments that a Fock space of 
flavor neutrinos does not exist. 

When neutrinos propagate through matter their behavior may be affected significantly, 
as was pointed out by Wolfenstein [12] in 1978 and emphasized by Mikheyev and Smirnov 
[13] as being important for neutrino oscillations. This is due to the fact that in the pres- 
ence of matter the effective mass induced by the forward scattering of neutrinos by the 
background changes the flavor oscillating parameters. A quantum mechanical treatment 
of neutrinos interacting with matter can be found in almost all books treating neutrinos. 
Peltoniemi and Sipilainen [14] have studied neutrino propagation in matter using a wave 
packet approach. Cardall and Chung [15] treats neutrino oscillations in a static uniform 
background by quantum field theory and show that they recover the quantum mechanical 
oscillation amplitude in the relativistic limit. Also Fuji et al. [11] treats neutrino oscillations 
in a static matter background using Bogoliubov transformations. 

In the present paper, we give a short presentation of the use of field theory methods to 
describe neutrino oscillations in vacuum and in matter. For neutrino oscillations in matter 
we carry out the calculations for the case of two flavors in detail, which displays the main 
features of our approach. Our aim is in particular to calculate explicitly the singularity 
structure of the neutrino Green's function in the presence of matter to study the resonace 
formula. 

In Sec. II, we review the basic formalism to be used. Section III contains an elaboration 
of the use of wave-packets for flavor states. We then calculate the oscillation amplitude for 
the case of vacuum oscillations in Sec. IV, and then, treat explicitly the case of two flavors 
in the presence of matter with constant density in Sec. V. 



The equation of motion for the neutrino field Uj, with mass rrij, can be written in the 
following form 



where we have introduced the source term xA x \ "which can be obtained from the complete 
Lagrangian density of the standard electroweak theory in the unitary gauge. The interaction 
part of primary interest to us is given by 



II. BASIC FORMALISM 
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where g is a dimensionless coupling constant and W a (x) and Z Q are the fields that describes 
the W and Z bosons with spin 1. The last part contains the Higgs field a and the electroweak 
vacuum expectation value v. We observe that the interaction part contains the auxiliary 
flavor fields vp. These fields are expressed in terms of the mass eigenfields uj by the relations 
vp(x) = Ylj UpjVj(x) and (pp{x) = ^ • UpjVj(x)), where U is unitary leptonic mixing matrix 
to be specified below. Originally the £7-matrix is defined only for the left handed fields. For 
simplicity we will assume that the same transformation also applies for the right handed 
fields. Inserting the expressions for vp and Vp expressed in terms of the mass eigenfields into 
the interaction Lagrangian, one obtains 
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We have omitted the part of the interaction Lagrangian containing the a field in the Higgs' 
mechanism. This Lagrangian gives 
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describing the fundamental processes that go on in the source and the detector. The first two 
terms in Eq. 3 define the flavor combinations through the intermediate boson interaction 
containing the U matrix components. This mixture in the interaction describes that the 
weak bosons create a coherent superposition of neutrino mass eigenfields at the vertex. At 
energy scales where the neutrino rest masses are small compared to their total energy, this 
can give rise to neturino oscillations. We also observe that the different neutrino fields 
Vj(x) are coupled to each other via the neutral current interaction. At the energy scale 
we are going to consider, the amplitudes describing the detection processes, e.g. in the 
SNO detector, only appear with the second order terms in the coupling constants. These 
processes comprise the elastic scattering of neutrinos on electrons, the inverse beta decay 
and the deuteron break up. To describe these processes we use the effective second order 
Lagrangian, obtained by neglecting the momentum of the intermediate bosons. This means 
that we insert the weak hadronic currents coupled to the neutrinos with the Fermi coupling 
constant Gp, and also use the analogue of this for the purely leptonic scattering in the 
Cherenkov detectors. Finally, we also use this effective Lagrangian to describe neutrino 
interactions with matter, neglecting recoil effects. The effective Lagrangian for neutrinos 
interacting with matter is 
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C(x)T tcI = ->/2G F £tf*^- L (^ (5) 

j=i V j =1 

where N e (x) and N n (x) are the electron and neutron density operators, respectively. In our 
further discussion, we will only consider the case when these densities are constant and take 
the neutrino fields to be the renormalized effective fields. 



III. THE FLAVOR STATES IN A WAVE PACKET APPROACH 

To zeroth order in the coupling constant g, the uncoupled Dirac equations for neutrinos 
(i$ — rrij)iSj(x) = 0, has the following free field solutions [16]: 



Vj\X) = 



Wf^/Wm £ («5(P)«i(P)^ X + bfbWWe**) , (6) 



where Uj and Vj are particle and antiparticle solutions to the Dirac equation. The conjugate 
field Vj(x) has an analogous expansion. This expression for the neutrino mass field in 
terms of the operators a|(p) and frf (p) can be inverted for the operators in terms of the 
fields in the usual way. In our view, only the mass eigenstate fields are those that build 
up the Fock space structure. The auxiliary "flavor fields" are defined as mixtures by the 
transformations given above, and are introduced for convenience. These flavor fields, of 
course, enter formally into the Lagrangian, as they represent the mixture produced by the 
weak interactions. However, they cannot be considered to represent asymptotic fields that 
are carrying irreducible representations of the Poincare or Lorentz groups. 

For our purposes we want to formulate the field theory slightly differently. We introduce 
the notation 

a' j (p,t)=a' j (p)e- iE >W i (7) 
where Ej(p) = \Jp 2 + rn 2 . The expansion can then be written 

" i{x) = I (2^3 v /2^(p7 ^ ( a ?( p ' *K(P) eiP " + *K(p)e- ip ") • ( 8 ) 
The time dependent operators af(p,t) can be projected out of the neutrino fields as 

af (p,t) = — L= iij(p) / d'xe^ujit^). (9) 

y/2hj(p) J 

These operators and their Hermitean conjugates create and annihilate particles of definite 
three-momentum p, at time t. 

A plane wave basis neutrino state \vj;r,p,t) with mass rrij, momentum p and spin 
r at time t is obtained by acting with a^(p,t) on the vacuum state |0), defined such 
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that (0|0) = 1. In order for \vj',r, p, t) to be covariantly normalized at equal times, i.e., 
(uj] s, q, t\uj] r, p, t) = 2Ej(p)(27r) 3 5^(p — q)5 rs , the definition for \vj\ r, p, t) is 

\u r ,r,p,t) = ^2E~{p)aJ(p,t)\0). (10) 

The auxiliary neutrino flavor fields z/ a , where a = e,fi,T, that we will use to construct 
our localized states, are expressed in the mass eigenfields Vj, where j = 1,2,3, by v a (x) = 
J2jU a jUj(x). The flavor fields are constructed with the mass eigenfields having the same 
three-momentum, p. The U^s are the elements of the unitary mixing matrix U, which for 
three generations is given in standard form by 

/ C 2 C 3 S 3 C 2 S 2 e~ i5 \ 

U = -S 3 Ci - S 1 S 2 C 3 e iS dC 3 - S 1 S 2 S 3 e iS S X C 2 (11) 
\ 5^3 - S 2 dC 3 e l& -S ± C 3 - S 2 S 3 C ie iS dC 2 J 

with Cj = cos 6i and Si = sin 9i, i — 1, 2, 3. When the CP- violating phase 5 = 0, the matrix 
elements are real and we have U*j = U a j for a — e, /i,t and j = 1, 2, 3. In what follows, we 
put 5 = 0. 

Let us define |f Q ;r, p, t) = - U*j\i>j; r, p, t)). Since we want to have an expression for 
\v a ; r, p, t) in terms of vj(x), we obtain 

k;r,p,*> = ^f/;. V / 2^(p)a; t (p,t)|0) = ^ W(p) / ^e^i/J^lO). (12) 

The state |i/ a ; r, p, t) has a definite momentum p, which means that Ap = 0. According 
to Heisenberg's uncertainty relation, ApAx > h/2, the state has therefore an infinitely wide 
spread in space. It also does not have a well-defined energy, since P° operating on the state 
does not satisfy E a = a/p 2 + m 2 a . This is of course one of the causes of the flavor oscillations 
in this mode of description. 

In order to have a Hilbert space state that is localized around the momentum P at 
position X and time T, one has to use a wave-packet with some distribution function 
F Q (X,P,p) = F a (P,p)e' iyi p . We can then write such a state \F a (r, P, X, T)), that is 
localized around P at the space time point X, T with spin projection r, as 

teP ' x ' r|>= /w A(P) K ' ,r ' p,n ( ' 

In what follows we will use the subscripts s and d to denote source and detector, respectively. 
Inserting our expression for \u a ; r, p, T), we finally have 

|F a (r,P„X„T,)) = / ^ ^^P) g^.„;(p) / ^P-x z/ t m)X)|0) (M) 

for £ = s,d. This expression can be taken as the definition of our localized neutrino flavor 
states. In actual applications, it is convenient to work with Gaussian wave packets, with a 
mean spread a around P. The corresponding wave function in the plane wave basis is 

/2tt\ 2 - (p ^ p)2 

F a (P e ,p)= l^j e £ = s ,d, (15) 
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where for simplicity will take the same ag for all flavors. The states in Eq. (14) have correct 
normalization for equal times if we make the interpretation E a (p) = U a jU*jEj(p), i.e., 
E a is the expectation value of the Hamiltonian in the state \u a ; r, p, T). 



IV. NEUTRINO OSCILLATIONS IN VACUUM 

The transition amplitude A r a s g = A r a s g(X. s: X d , T s , T d , P s , P d ) in the rest frame of the 
detector and the source for a neutrino of flavor a, momentum around P s , and spin r to be 
produced at the source and to be detected as a flavor f3, with momentum around P d , and 
spin s is given by 

A r a s p = {Fp(s,-p d ,X d ,T d )\F a (r,-p s ,X s ,T s )). (16) 
Using the expression for \F a (r, Pg, Xg, Tg)) and its conjugate, we have 

x(0\u j (T d ,y)vl(T.,x)\0). (17) 

The expression (0\vj(y)vl(x)\0) is the propagation amplitude for a mass eigenstate k being 
created at x — (T s ,x), propagating to y — (T d ,y), where a mass eigenstate j is being 
annihilated. In vacuum, j has to be equal to k. On the other hand, when neutrinos 
propagate through matter it is possible that j ^ k. Since T d > T s i.e. y° > x°, we can insert 
the time-ordering operator T obtaining 

(0\u,(y)ut(x)\0) -> (0\T(u,(y)^))m 3k . (18) 

The last expression is simply 

(OFteMvfam = (0|T(^(y)F J (x))|0) 7 ° = S 3 (y - x)^, (19) 

where Sj denotes the Feynman fermion propagator defined by 

^-*) = '/S^Tfc e ^ > - (20) 

The integration in the complex A;°-plane has to be taken along the whole real axis. We can 
then rewrite the expression for the amplitude as follows 

ats = vr; rr ! d * q I d " p F *^ p ^^ F a (x s ,p s , P ) 

x J d 3 ye~^ J d 3 xe^* J _0_ e -^(T d -T s )-k.( y -x)] 

x . uf(q)^ + m j )-f°u r j (p) 
k 2 — m 2 + ie 
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We next calculate the integrals in the middle in the order x, y, and k. This gives 



3 ip-x (f + m j) r -ik(y-x) 
7 9 2 i • 



(22) 

Closing the contour in the lower k° half-plane and calculating the residues of the positive 
energy poles, it follows that 

ats = vr/ rr I d " q I d ' p F ^ p i^) F a (x s ,p s , P ) 

a " V * * 7 (2-) 3 J (2vr)3 ^/s^J ,/2K(p7 

x(2tt)V 3 )( P - q ) e -^-W(p) st (g) + ^l /^(p). (23) 



We finally obtain 

d 3 p F*(X d ,P d ,p)F a (X s ,P s ,p) 

3 



2tt) 3 v^pj v^(pT 

x <(p) jjj + m 3 -) 7 % r (p) e _^. (p)(rd - Ts) _ (24) 
2Ej(p) 

The part (p) (^ + rrij) j°Uj(p) can be simplified by using that (J> — m,j)Uj(p) = 0. This 
gives 

uf(p) (j> + m 3 ) 7°«;(p) = (2^(p)) 2 ^. (25) 
We finally arrive at the following rather simple expression for the transition amplitude 

A "" - J w? jmm (26) 

where 

= 5" E ^^2^(p)e-^^^-^). (27) 

3 

The amplitude in Eq. (26) will be our standard reference amplitude, when we later 
calculate the amplitude for neutrino oscillations in the presence of matter. Using Gaussian 
wave packets, it is possible to show that we regain the usual neutrino oscillation formula 
in the ultra-relativistic limit. In particular, it is useful to consider the time average of the 
amplitude, as we do not know when the neutrinos are emitted. We will not here elaborate 
further on this point, which is treated in Ref. [10], but continue with the case when neutrinos 
interact with matter. 



7 



V. NEUTRINOS IN MATTER: THE RESONANCE FORMULA 



A. General discussion 

We will next proceed to derive the neutrino oscillation probability for neutrinos prop- 
agating in matter with constant density. Although a perturbative treatment of neutrino 
oscillations in matter would give the influence of matter for density factors that are small 
compared to the energy differences Am 2 /E between the neutrinos, the behavior of the neu- 
trinos are quite different for densities that are of the same order as Am 2 /E. This case is 
known as the MSW effect and was pointed out by Wolfenstein, Mikehyev, and Smirnov 
as being potentially important. The quantum mechanical case can be found in many text 
books. To treat this effect with quantum field theory we have to calculate the Green's 
function for the neutrinos in matter. The equation of motion derived from C™ ttcr (x) is in 
component form 

- mj )v 3 {x) + J2 4jMx) = Xj (x), (28) 

k 

where Xj( x ) i s the source. Here the term responsible for the interaction with matter is given 
by contributions from neutral (N) and charge changing (C) currents in the static limit as 

A»(x) = ^(AnI + A C K), (29) 

where Kjj, = U*jU e k is the projector onto the electron neutrino flavor space. The terms A^ 
and A c have the values A N = — y/2GpN n /2 : and Ac = y/2GpN e expressed in terms of the 
number density of neutrons or electrons, respectively. For antineutrinos there is a change of 
sign in these terms. 

By transforming to momentum space and assuming that the neutrinos all have the same 
three-momentum p, the Dirac equation in the presence of matter can be written in momen- 
tum space as 

(#7° — p • 7I — M — (A N 1 + A c K)^°)u(p) = 0. (30) 

In this section, we will use the symbols um for the solutions to the Dirac equation in the 
presence of matter. 

The Eq. (30) can be written as 

Hu M {p) = (-P • 77° + M7 + A N t + A c K)u M {p) (31) 

To obtain the full set of roots one has, in the energy representation, to enlarge the flavor 
space with the two-dimensional energy space, by using u± = |(1 ± 7°)w for each flavor /. 
This gives a 2f— dimensional representation of the Hamiltonian, that will have 2f roots, i.e. 
/ roots of both signs. These roots are calculated for / = 2 later. 

It is then possible to diagonalize the Hamiltonian H, which is symmetric in this 
2/— dimensional representation, with a unitary transformation Um which is analogous to 
that one given e.g. in Ref. [17] for the case of an ordinary quantum mechanical description. 
The Eq. (31) then becomes 
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Hu M j = EjiiMj, (32) 

where um is related to u in the mass representation by % = UmU^u. Since the flavor 
states Uf are related to the matter states by Uf — U m um, where U is naturally extended 
to the same dimension, the flavor states can be expanded in the mass representation also 
in case of neutrinos propagating through matter (of constant density). For the antineutri- 
nos a corresponding equation holds with different sign on the matter terms, and the same 
transformation Um will diagonalize the Hamiltonian also for the antineutrinos. 

The Feynman propagator, as a Green's function, is invariant under unitary transforma- 
tions, and can therefore alternatively be expressed in the mass representation. 

To obtain the transition amplitude we now utilize the same approach as in the previ- 
ous section. The neutrino field is expanded in the solutions to Eq. (30) and the spatial 
momentum plane waves. The corresponding annihilation operator satisfies the equation 

a a j (p,t)=a s j (p)e- iE i t , (33) 

where Ej now is the positive eigenvalue of H and the two negative eigenvalues go with the 
adjoint creation operators. For the antiparticle operators these change their roles and the 
sign of the matter terms also change. 

The amplitude of interest is in this case given by the expression 

A rs f d*q f d*p F;(X d ,P d ,q)F a (X s ,P s ,p) 

A °* J (2k)*J (27T)3 ^2EM y/^M ^ 

where 

MH aP (q,p) = J2 U MP 3 U Mak M^ jk (q,p) 

jk 

= £ U M(3 jU Mak J d 3 x J r/V" x ' " i y 1^,(11 - x), (35) 



jk 

and the expression R M jk( x ) i s gi ven by 



R r Mjk{*) = J ^e-^(g) ?7 ^0|T(z,,(y)^( a ;))|0) 7 V Mfc (p). (36) 

Using the transformations given above and the fact that the time-ordered Green's function 
can be expressed in any basis, and thus can be replaced by the Feynman propagator in the 
mass representation, we arrive at the expression 

M r ^(q, p) = £ U Pj U* ak J d 3 x J tfy ip - x e- iqy i^(y - *), (37) 

jk 

where 

Rji(x) = J 7^I e "^(^)^°^(^,A;) 7 V( P ). (38) 
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The expression M^ a/3 in Eq. (37) will be our main task to calculate. 

We first want to calculate the Feynman propagator in the presence of matter. What we 
seek is thus the inverse S(A, k) of the operator 

S~ 1 (A, k) = G7 — k • 7I — M, (39) 

where G = k°l — A N 1 — A C K which operates in the tensor product space of the four 
dimensional Dirac space with the three dimensional mass space in the mass representation. 

The neutral currents give a constant contribution to the energy that can be absorbed in 
its definition until the end of the calculation and will therefore be neglected from now on. 
The charge current contribution, on the other hand, gives a more complicated contribution 
and will be dealt with explicitly. For convenience we will use the notation A = A c below. 
The calculation can then readily be done and one finds 

S(A, k) = [G7 - k • 7I + M + (G7 - k • 7I + M)^N^C] D~\ (40) 

where N = G 2 - k 2 t - M 2 , C — -A [M, K] , and D = N - CN~ l C. This means that we 
are interested in the poles of D -1 with respect to k°. We will perform the analysis for the 
case of two flavors below, leaving the three-flavor case for the future. 



B. The case of two flavors 

Here we will confine ourselves to calculate the amplitude explicitly only for two flavors, 
since we can then easily compare our result with the quantum mechanical expression for the 
resonance formula. 

For two flavors we thus want to calculate the appearance amplitude 



d 3 g r d 3 p F;(X d ,P,,q)F e (X 5 ,P s ,p ) 
(2*)* J (2tt)3 v /2^(q7 y^EJv) 



^matter - / ,r> x 3 / x 3 , t-t , . r M Me^V)- {^) 



Carrying out the integrations as before we end up with the expression 
A rs _, d 3 p F;(X d ,P d ,p)F e (X s ,P s ,p) 



dk° 



x i I ^ e - tt0 C r --^«'(p)7°5(A {k°, P)) 7 °<(P). (42) 
We are thus interested in the expression 

MZe, = * J ^e-^u^SiA, (k°, p)) 7 °<(p), (43) 

where ue = Y^=i ^tj u 3i ^ = e, and T = Td — T s . 

In the case of vacuum propagation, we have M™ = 5 rs ^ U*jUnj2Ej(p)e~ lEjT , which in 
the high energy limit goes to 5 rs e~ iET 2E sin 26 sin[T(£ 2 - E 1 )/2], with 2E = E 1 + E 2 . In 
what follows, we will frequently omit the spin indices. 
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Let us first calulate the matrix D = N — CN 1 C for the two flavors e and The 
calculation is done in the mass eigenstate basis. In this case, we have 

G = k°l 2 - AK. (44) 

Remember here that k° is really k° — A^. Then, 



C = -A [M, K] = A(m 2 - m,) ( K ™ ) 



(45) 



where K 21 = K 12 = sin 9 cos 9. Thus, 



c = K {-io) («) 

with k = A{m 2 — mi)Ki 2 . We can then readily calculate the expression CN~ l C . The result 
is 

CN-'C = —N T = —N, (47) 

detAT detAT ' y ! 

since A" is symmetric. Thus, 



k 2 , T ( ' 



detN I detiV 1 



and 



D-i = - ( ^ ^ ] (49) 



In calculating the expression wy°S(A, (k°, p))7°-u, we can use the fact that the spinors u 
fulfill the relations 

u(p)((E- AFsr) 7 ° - p • 7 ]1 - M) = 0. (50) 
This leads to two different terms 

S (1) = (k° + E - 2AK)D- 1 (51) 

and 

S (2) = (k° + E- 2AK)N- 1 CD- 1 (52) 

with 

u a °S(A, (k°, p))^ ^ = u\sV Uj + tiiSffuj. (53) 

We are looking for the singularities of D~ l and N~ 1 CD~ 1 . The singularities of D~ l are 
given by the roots with respect to k° of the equation 

N n N 22 - Nl 2 + k 2 = 0. (54) 

Now, A" can be calculated to be 

M _ ( (k ) 2 - p 2 - m\ + (A 2 - 2k°A)K n (A 2 - 2k°A)K 12 \ 

iV _ I (A 2 - 2k°A)K 12 (k ) 2 - p 2 - m 2 + (A 2 - 2A;°A)K 22 J ' 
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The equation (54) for k° is then 

((A; ) 2 - p 2 - m\ + (A 2 - 2k°A)K u )((k ) 2 - p 2 - m\ + (A 2 - 2k°A)K 22 ) - 

(A 2 - 2k°A) 2 K 2 12 + k 2 = . (56) 

The other singularities are determined by the expression N^CD^ 1 . This can easily be 
calculated to be 

N^CD- 1 = -— £ ( ° I V (57) 

Obviously, the singularities are again at the same location as the roots of Eq. (56). 

To solve Eq. (56) we first simplify the equation using the properties K\\K 22 — K\ 2 = 
and Ku + K 22 = 1, since K is a projector. We also use E 2 = p 2 + m 2 for % = 1,2. Then 

((A; ) 2 - E 2 )((k ) 2 - E\) + A(A - 2A;°)((A; ) 2 - K XX E\ - K 22 E 2 ) + k 2 = (58) 

Finally, we also have k 2 = A 2 (Am) 2 K u K 22 = A 2 (Am) 2 \ sin 2 26. The term K n E% + 
K 22 E\ can be simplified to be \{E\ + E 2 ) + \(m\ — m\)cos26. If we, for typographical 
reasons, put k° = x the equation can be written 

x 4 - 2Ax 3 + (A 2 - El - El)x 2 + A(El + El + (m 2 - m\) cos26)x 
-\a 2 (E 2 x + El + (m 2 2 - ml) cos 26) + ^A 2 (Am) 2 sin 2 26 

+ElEl = 0. (59) 

We further introduce new variables E = (E\ + E 2 )/2 and Am 2 = m 2 , — m\. Then, the 
equation above can be written 

x * _ 2Ax 6 + \A z -2E z -2[ ] ] x 2 + A | 2E 2 + 2 ( + Am 2 cos 20 j x 

+ A//,* cos 20 \ A 2 ( A '"~ ) sin 2 26 

\2(mi+m 2 )J 

, 2 (Am 2 \ 2 ( Am 2 V , s 

+e ' 2E {^e) + {^e) = ' ^ 

which can be solved exactly. Let us call the roots of this equation Aj, i — 1, ... ,4. 

These roots are of course exactly the same as those we get by solving the secular equation 
for the Hamiltonian given in Eq. (31), when we enlarge the two-dimensional flavor space with 
the two-dimensional u± space related to diagonalizing 7 , as was mentioned in Sec. VA. 

Since we want to calculate the amplitude MJ* we are, in principle, interested in the 
expression 

U^S^Ue = -SCUi7° 1 S'i 1 7 -U 1 - S 2 Mi7°S'i27 -U2 + ^U^S^^Ux + CSU 2 ^f° S 22 ^f°U 2 

= sc(u 2 ^°S 22 ^°u 2 - ■ui7°S'ii7 -ui) + (c 2 -u 2 7°»S'2i7 Mi - s 2 ui^° Si^u^, (61) 
where s = sin0 and c = cos 6. 
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The amplitude can be written in the form 



M£ = = S rs det ^ +K2 ((k°) 3 a 3 + (kya 2 + *% + a ), (62) 

where the coefficients a, have to be calculated. We will discuss this point later and first 
calculate the amplitude. 



C. Calculation of the amplitude 

The roots \, % — 1, . . . , 4 of Eq. (60) are all real, and we will proceed with the calculation 
of the amplitude. We will also assume that the roots are indexed so that in the limit A — > 
the roots take on the values \± = E±, \ 2 = E 2 , A 3 = —Ei, and A4 = —E 2 . We can then 
close the contour in the lower half-plane and calculate the residues of the amplitude avoiding 
the negative roots. The amplitude is then (minus) the sum of the residues of the poles at 
Ai and A 2 . 

The sum of the residues of the amplitude M eil is then 

y, ■ A^a 3 + A| a 2 + AjOj + Qp _iTAi _ ■ A 2 a 3 + A 2 a 2 + A 2 ai + Qq -«ta 2 / fi o\ 

e " = '(A 1 -A 2 )(A 1 -A 3 )(A 1 -A 4 ) e ~\A 2 -A 1 )(A 2 -A 3 )(A 2 -A 4 ) e ' [ } 

This expression is similar to that of the matrix element of the spectral resolution of the 
evolution operator, only that here there is an influence from the negative energy terms 
through the denominators. 

We first introduce the difference of the roots as A = A 2 — Ai. Next, we observe that 
we can translate all the roots by A/2 so that x = y + A/2 to make the sum of the new 
roots yi,i — 1, . . . , 4 satisfy y\ + y 2 + y% + 2/4 = 0. This transforms away the x 3 -term in the 
equation. Also the denominators in the amplitude are unaffected by this translation. The 
numerators will obtain terms proportional to A or A 2 , which we can safely neglect. 

In calculating the coefficients a,, % — 0, . . . , 3 we will neglect small terms in the numera- 
tor. Our main interest is to locate the resonance region where the denominators are small. 
We therefore first put the matter parameter A to zero in the wave functions as well as in 
the terms S^ l \i = 1,2. We further use the approximations ufu r j S rs (Ei + Ej) for all 
and u^vJj ps 5 rs (mi + mj). For the diagonal elements the equality is strict. 

In the limit of neglecting small terms in the numerator, we obtain the compact formulae: 

a 3 ps sin 29 Am 2 — , (64) 
2E 

a 2 sin 26* Am 2 , (65) 
E 

a 1 ^sm29Am 2 — , (66) 

and 

a = 0. (67) 
The large contributions to the a^s given above all come from the term. 
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We now realize that introducing the small shift in the energy by the terms An, as we did 
in the beginning, amounts only to add An to A in the roots of Eq. (60). This is still a very 
small term and does not change anything essential in the above approximation, especially 
not in the denominators, that contain differences of the roots. Apart from appearing in the 
small terms neglected, the contribution from An appears in the overall phase factor as an 
addition to A. However, this will not affect the probability. 

Finally, we approximate the sums and products of roots in terms of yi with their limiting 
values as A = (and A N = 0). This leads to an expression for the amplitude of the form 

Me " ~ -i i ^e-' T <«' + »- 4 >/ 2 j-i (a 3 2B 2 + 2Ea 2 + 2a, + ^) sinlTA/2) 

+ f^ + !^-!^W4/2)l. (68) 



2E 2E 3 

Inserting the values for the coefficients from Eqs. (64) — (67), we obtain the formula 

~ -t5 rs e~ i{E - A/2)T 2E sin 26 M cos(TA/2) - i sin(TA/2)^ , (69) 

where 

. nn Am 2 sin 29 

sm26 M = . (70) 

In the limit A = 0, we have A = Am 2 /2E. Neglecting the term of order Am 2 /E 2 , we obtain 

M e/I (A = 0) = e~ iET 2E sin 26 sin(TAm 2 /4£), (71) 

which is the previously derived result without matter effects and coincides with the quantum 
mechanical expression. 

We can simplify this result in Eq. (69) further by observing that A = {y\ — y\)j2E within 
the approximation scheme we are envisaging. We will use a perturbative approach to obtain 
approximate analytic expressions for the roots y^. First, we omit the linear term in Eq. (60) 
after translating to x = y + A/2 and solve the resulting equation in y exactly. We then use 
dimensional analysis to make the Ansatz 

yi (e) = yj ^A 2 + B 2 - VA 2 B 2 - A 2 S 2 + 4a 2 - AAm 2 cos 26e, (72) 

y 2 (e) = yj ^A 2 + B 2 + VA 2 B 2 - A 2 S 2 + Aa 2 - AAm 2 cos 26e, (73) 

where B 2 = E 2 + {^r) 2 and a = Am 2 /4. Inserting this into Eq. (60) gives the equations 

e = yi {e), i = 1,2. (74) 



Solving to lowest order gives e\ = \j B 2 — 2a = E 1 and e 2 = \/ B 2 + 2a = E 2 . The difference 
of the roots is then 
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A m -^-^ (^\Ja 2 B 2 - A 2 S 2 + 4a 2 - AAm 2 cos 26V B 2 + 2a 
+ \/ A 2 B 2 - A 2 S 2 + 4a 2 - AAm 2 cos 26 4 B 2 - 2a^ . 



This gives 



where 



sin 26 M 



Am 2 sin 26 2 sin 26 



2EA 



A 



(75) 



(76) 



A = \J {AB/2a) 2 - {AS /2a) 2 + 1 - (A/ a) cos26VB 2 + 2a 
+ \J (AB/2a) 2 - (AS/2a) 2 + 1 - {A/a) cos26V~B 2 - 2a. 

The minimum of A is given up to small terms by 

cos 26=— (VB 2 - 2a + VB 2 + 2a) = ' " ' 



4a 



Am 2 



(77) 



(78) 



This minimum location coincides with the quantum mechanical location. However, the 
minimum value is not exactly the same as in the quantum mechanical case due to small 
extra terms. One of these terms is dependent upon the sum of the neutrino masses rather 
than their mass squared differences. 

The expression for A can be rewritten as 

A = \J (AVB 2 + 2a/2a - cos 2#) 2 + sin 2 26(1 - A 2 /( mi + m 2 ) 2 ) - A 2 /2a 

+ \J (A4W^2~a/2a - cos2fl) 2 + sin 2 26{\ - A 2 /{m x + m 2 ) 2 ) + A 2 /2a. (79) 

This is the resonance formula. The minimum at the resonance peak is dependent on the 
neutrino mass in the very small terms: 



A pcak = 2 sin 20 ( 1 



A 2 



2 (mi + m 2 y 



(80) 



With realistic values of A pz 10~ 12 eV, a pz 10~ 5 eV 2 , and B in the MeV range or higher, 
consistent with the experimental mass differences and the possible matter densities, the 
corrections are of the order of 10 -14 , and the value of A at the minimum is given by A = 
2v/sin 2 20(1 + 0(A 2 /(m x + m 2 ) 2 ) + 0(A 2 / Am 2 ) pz 2 sin 26 as in the quantum mechanical 
case. This indicates that the other corrections would be still smaller. The full expression 
for A can be calculated using the formulas above without approximations, but will be quite 
difficult to handle and is not of much use. Our expression Eq. (69) contains the field theoretic 
roots to the equation. These roots are dependent not only on the mass squared differences, 
but also weakly dependent on the sum of the neutrino masses. This comes from the term k 
in the equation for the singularities. This term appears from the nontrivial coupling of the 
two mass eigenstates in the calculation of the inverse propagator. 
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VI. SUMMARY AND CONCLUSIONS 



In summary, we have discussed field theory methods to treat neutrino oscillation phe- 
nomena with wave packets for neutrino flavor states in vacuum and in matter with constant 
density. We have carried out the calculation of the oscillation amplitude with matter effects 
for the case of two flavors explicitly. The resulting formula coincides with the quantum 
mechanical expression using plane waves, up to very small terms. These terms depend on 
the neutrino masses, but are unfortunately too small to be observable at the present time. 
To the best of our knowledge this is the first detailed calculation of the neutrino oscilla- 
tion formula using field theory and integrating out the poles in the Green's function of the 
neutrinos in the presence of matter. Earlier works, e.g. in Ref. [11,15] have used various 
approximations at an earlier stage and thereby not being able to isolate the poles exactly. 

Our calculation is done in what is essentially the rest frame of the detector and source, 
using wave packets, in which the neutrinos have the same three-momenta. It should be 
possible, though more complicated, to relax this condition, and let also the three-momenta 
in the superposition vary within the limits allowed by the wave packet parameters. 

For three flavors the corresponding equation for the poles will be of sixth order, and 
can therefore not in general be solved exactly. We nevertheless believe that approximate 
analytic solutions to the equation can be found that can be used to study the singularity 
structure, in a way analogous to the two flavor case treated here. 
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